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ABSTRACT 


As part of the Jet Propulsion Laboratory's Low-cost Solar Array 
Project, an analytical method has been developed for determining the 
minimum thickness for simply supported, rectangular glass plates 
subjected to uniform normal pressure environmental loads such as wind, 
earthquake, snow, and deadweight. The method consists of comparing an 
analytical prediction of the stress in the glass panel to a glass 
breakage stress determined from fracture mechanics considerations. 
Based on extensive analysis using the nonlinear finite element 
structural analysis program ARGUS, design curves for the structural 
analysis of simply supported rectangular plates have been developed. 
These curves yield the center deflection, center stress and corner 
stress as a function of a dimensionless parameter describing the load 
intensity. Results are included for plates having length-to-wid th 
ratios of 1, 1.5, 2, 3 and 4. The load range considered extends to 
1000 times the load at which the behavior of the plate becomes 
significantly nonlinear. Over the load range analyzed, the analysis 
shows that the ratio of center deflection to plate thickness for a 
plate of length-to-width ratio of 4 is less than 70 to 1, whereas 
linear theory would predict a center deflection about 1200 times the 
plate thickness. The stress is also markedly lower than would be 
predicted by linear theory. These analytical results show good 
agreement with the analytical and experimental work of others. 

A method of estimating the glass breakage stress as a function 
of a specified failure rate, degree of glass temper, design life, load 
duration time, and panel size is presented. Development of this 
method consisted largely of collecting and/or adapting, in convenient 
form, the best available information from the literature. To 
establish the glass breakage stress versus probability of failure, the 
experimental data of other investigators has been reanalyzed to obtain 
a "best-fit" Weibull statistical distribution. This state-of-the-art 
analysis yields the glass breakage strength as a function of failure 
probability. 


A step-by-step procedure is given which includes determination 
of the applied stresses, glass breakage stress, and criteria for 
comparing the applied stresses to the estimated glass breakage 
stress. Errors in the method for determining the applied stress 
should be largely self-correcting, since this same method has been 
used to convert the glass plate burst pressure data of various 
investigators to the glass breakage stress reported herein. The 
results of this analysis indicate whether or not an assumed glass 
thickness will withstand the design loads at a failure rate specified 
by the designer. Finally, a sample problem is presented in Appendix I. 
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SECTION I 
INTRODUCTION 


The need for a straightforward method of determining the glass 
thickness icquired for solar collector panels is evident. These 
panels are subjected to environmental loads which may be characterized 
as uniform normal loads (wind, earthquake, ice, snow and deadweight). 
Currently, the sizing of rectangular glass windows for uniform 
pressure load 9 is based largely on design curves available from glass 
manufacturers. These curves show the glass thickness versus window 
area required to sustain a given windload with a failure rate of 8 
windows per 1000. These design curves are based on an empirical 
equation of the form 

P = k (t+t 2 ) 

where P is the total load, t is the thickness, and k is an empirical 
constant. Empirical equations of this type have been developed by Orr 
(Reference 1) and others. The application of this method for 
determining the thickness required for glass in solar collector panels 
has several shortcomings: 

(1) Since relatively few samples of any particular size and 
thickness were tested (Reference 1), the resulting design 
curves do not have a strong statistical base. 

Statistically significant sample sizes are known to be 
important when dealing with brittle materials, such as 
glass . 

(2) Loading was applied incrementally (Reference 1), and each 
load increment was held constant between 5 and 25 minutes 
while deflection, curvature and strain measurements were 
taken. Consequently, the load-time history at fracture 
for the various panels is not known. As will be seen 
later, the breakage strength of the glass depends on the 
load duration time. 

(3) Because solar collector panels are tilted with respect to 
the horizon, loads such as snow and deadweight are 
applicable to the design of photovoltaic panels. Such 
loads, however, are not necessarily important factors in 
the design of ordinary windows. Because these loads are 
of longer time duration than wind loads, they will cause 
failure at a lesser applied stress, and therefore require 
special consideration in determining the required glass 
thickness . 

It should be pointed out that since the above reservations 
concerning current methods for sizing rectangular windows were 
enumerated, the Pittsburg Plate Glass Company has published an 
improved glass thickness selection method entitled Glass Thickness 
Recommendations to Meet Architects Specified 1-Minute Wind Load" 
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(April 23, 1979). PPG's new method is soundly based on the same 
fundamental method of structural analysis that is used in this 
. report: comparing a calculated applied stress to the allowable 

| strength for glass. PPG also uses nonlinear finite element stress 

analysis to obtain an accurate estimate of the applied stress. 
However, the scope of the above document is. limited; their glass 
thickness recommendations are for wind loads of 1-minute duration and 
a failure rate of 8 windows per 1000 at the design load. 



SECTION II 


FAILURE OF GLASS 


Weiderhorn (Reference 2) has estimated that the theoretical 
cohesive strength of silica glass is of the order of 3 x 10^ psi. 

The highest measured values for the strength of glass (glass fibers 
tested in a vacuum) approach this value. On the other hand, values 
normally associated with the engineering or practical strength of 
glass are approximately 10 to 1000 times lower than the theoretical 
cohesive strength. This discrepancy between the theoretical and the 
practical strength of glass is attributable to the brittle nature of 
the material. Stress concentrations at the tips of existing flaws, 
especially surface flaws, are not relieved by plastic yielding of the 
material. For this reason glass always fails in tension when these 
flaws grow (under sustained loading) to some critical size. Thus, the 
surface condition of the glass is the most important single factor 
influencing the breakage strength of any glass part. The following 
factors are known to materially influence the practical strength of 
glass parts. 

1. Part Size . Even though the intrinsic strength of glass does not 
vary with part size, small glass parts (glass fibers, for example) 
exhibit strengths very much higher than relatively larger parts. 

Large parts have more and larger flaws, increasing the probability 
that a severe flaw will coincide with a region of high tensile 
stress. This phenomenon leads to a lower breakage strength for larger 
glass parts. 

2. Moisture . Glass parts which have been dried and tested in a 
vacuum exhibit higher strengths than those tested in the presence of 
moisture. This behavior has been predicted from the "stress corrosion 
theory.” It is important to note that the small amount of moisture 
normally present in the atmosphere is sufficient to account for most 
of this effect. Therefore, it will not be necessary to differentiate 
between the strength of glass parts in wet versus dry atmospheric 
conditions during normal outdoor exposure. 

3. Surface Damage . The strength of glass parts may be expected to 
deteriorate with time due to environmentally induced damage resulting 
from handling, shipping, installation, cleaning, hail and rock impact, 
sandstorms and other causes. Dalgliesch and Beason (References 3 and 
4) have indicated that there is evidence of a reverse effect as well 
in which severe stress concentrations at the tips of sharp microcracks 
are supposedly alleviated by rounding off those cracks for which the 
stress corrosion is not made highly directional by a large applied 
stress . 


4. Du ration of Load . Glass parts will sustain loads for short 
per iodT”of time which will result in eventual failure. The existing 
flaws in the glass grow under the influence of a sustained load. When 
these cracks reach some critical size, failure occurs. 


2-1 


SECTION III 


OVERVIEW OF PROPOSED METHOD 


The shortcomings of existing methods of determining the required 
thickness for simply supported, rectangular glass panels subjected to 
uniform normal pressure loads have already been noted. The existing 
methods are based on experimental burst pressure data for glass of 
various sizes and thicknesses, thus circumventing the need to predict 
the applied stress. The essence of the proposed method follows the 
conventional theme of structural analysis: applied stress is pre- 

dicted as accurately as possible and subsequently compared to an 
allowable stress. The allowable strength will be based on those 
factors known to influence the breakage strength of glass: degree of 

temper, duration of the load, and size of the glass part; Addition- 
ally, the strength of the glass is treated in a probabilistic sense, 
so that the allowable strength is actually the glass breakage strength 
at a specified probability of failure. The proposed method is 
depicted graphically in Figure 1. 



Figure 1. Diagrammatic Representation of Proposed 
Glass Thickness Sizing Method 
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In order to Implement the proposed glass thickness sizing method 
the following steps were taken: 

(a) A nonlinear finite element analysis of simply supported, 
rectangular plater? subjected to uniform normal pressure 
loads was performed over a wide load range for plates with 
length-to-width ratios of 1, 1.5, 2, 3 and 4. The results 
of this analysis are presented in nondiraens ional nomograph 
form. 

(b) Means to quantify the effect on glass strength of degree 
of temper, duration of loading, and size of glass plates 
have been adapted from the literature, 

(c) The nondimensional design curves are used to convert the 
glass pressure burst data of various investigators to 
glass breakage stress. This procedure minimizes errors 
resulting from inaccuracies in the proposed method for 
estimating applied stress, since consistent errors exist 
in the glass breakage stress reported herein. 

(d) Weibull statistical analysis has been performed on the 
extensive test data of Beason (Reference 5) and Bowles and 
Sugarraan (Reference 6). These investigators performed 
burst pressure tests on annealed sheet and plate glass of 
various thicknesses and sizes. The results of this 
analysis give the strength of annealed glass versus 
expected failure rate. 

(e) Additionally, Abraham Wilson (Reference 7) of the Jet 
Propulsion Laboratory has performed burst pressure tests 
on 48 x 48 x 0.125 inch simply supported, tempered glass 
plates. The results of a Weibull statistical analysis of 
these data give the strength of 0.125 in. -thick tempered 
glass versus the expected failure rate. 

The development of this methodology is discussed in the following 
sections . 
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SECTION IV 


CAVEAT o 


The proposed method is intended to aid the practicing 
professional engineer or architect in determining the glass thickness 
required for a solar collector panel for a given application. The 
method is based on state-of-the-art analytical techniques (nonlinear 
finite element analysis and Weibull statistical analysis). While the 
method is believed to be sound, the author cannot assume 
responsibility for any loss incurred as the result of failure of a 
glass part designed using this method. Limitations of the proposed 
method recognized at this writing are enumerated below. 

(1) It is limited to rectangular glass panels which are simply 
supported (negligible inplane or rotational constraint) on 
all four edges. This limitation does not apply if the 
stress in the glass panel can be assessed accurately by 
other means. 

(2) It is limited to those loading conditions which can be 
considered to be uniform normal loads, such as wind, 
earthquake, snow, and deadweight. 

(3) It is not intended to cover localized impact load', such 
as those resulting from hail or other missiles, lor 
guidance in designing hail-resistant glass solar collector 
panels, the reader is referred to "Photovoltaic Solar 
Panel Resistance to Simulated Hail" (Reference 8). 

(4) It does not cover thermally induced stresses. These might 
arise from differences in the thermal expansion 
coefficients of glass and its mounting frame, from uneven 
temperatures across the surface of the panel due to 
shading from the adjacent row(s) of panels, or from local 
temperature gradients attributable to back-biasing of 
damaged solar cells. 

(5) It does not account analytically for deterioration of the 
strength of the glass with time resulting from handling, 
shipping, installation, cleaning, missile impact, and 
other environmentally induced damage. These effects 
cannot be assessed quantitatively at this time. On the 
other hand, the glass breakage data analyzed in Section 
VI, which include 20-year-old weathered glass samples, 
show that while the mean strength of weathered glass is 
less than that of new glass, the breakage strength for a 
1% probability of failure is about the same. 
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SECTION V 
STRESS ANALYSIS 


A good estimate of the stress in a uniformly loaded, simply 
supported glass plate is a necessary ingredient of the proposed method 
to determine the required thickness of the glass. It is well known 
that plates loaded normal to their plane develop substantial mid-plane 
membrane stresses at moderate loads. This results in a nonlinear 
stiffening effect as the plate is loaded beyond the point at which the 
center deflection is about one-half the plate thickness. The stress 
and deflection of plates laterally loaded at moderate to high loads are 
significantly lower than are predicted by the linear theory, which 
assumes that the lateral load is resisted by bending stresses only and 
does not include the stiffening effect of the membrane stresses 
developed as the plate deflects under load. In practice, a linear 
analysis of glass plates subjcted to normal loads results in an 
over ly-conservat ive design; that is, the glass is thicker than 
necessary . 

Existing analytical and experimental results reviewed by the 
author (References 6, 9, 10, 11 and 12) have one or more of the 
following limitations with regard to their usefulness to the designer: 

(1) They refer to glass whose edges are not simply supported 
( i . e . , free to rotate but not to slide in plane) 

(2) Only one length/width ratio is considered 

(3) They are not presented in dimensionless form 

(4) Only a narrow load r^nge is included 

Therefore, a nonlinear stress analysis of simply supported, rectangular 
plates subjected to uniform normal pressure loads is required. 

A. NONLINEAR FINITE ELEMENT z\NALYSIS 

To obtain a good estimate of the stress in simply supported, 
rectangular glass plates subjected to uniform normal pressure loads, an 
extensive nonlinear analysis was carried out using the ARGUS nonlinear 
structural analysis program (Reference 13). The following symbols are 
used in reporting the results of this analysis and throughout the 
report in calculating the stresses in a simply supported, rectangular 
plate subjected to a uniform normal pressure load. 

a = length of plate 

b = width of plate 

t = thickness of plate 

E - Young's modulus of plate material 
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v = Poisson's ratio of plate material 
D = Flexural rigidity of plate 



12(l- V 2 ) 


“ center deflection of plate 

o = maximum positive principal stress in plate 

The maximum positive principal stress (o) is taken as the 
critical stress, since glass is a brittle material. Figure 2 shows a 
simply supported, rectangular plate subjected to a uniform normal 
pressure load and defines the coordinate system used. 

Plates with length/width ratios of 1, 1.5, 2, 3 and h have been 
analyzed. This analysis has been carried out over a broad loading 
range, extending to 1000 times the loading for which the results 
become significantly nonlinear. The finite element models of plates 
analyzed showing the various length/width ratios are shown in Figure 
3. Note that due to symmetry only the upper right quadrant of the 
plates is modeled. 

In general, the more finely divided a finite element model is, 
the more accurate the results. However, making the grid finer 
increases analysis costs geometrically. Nonlinear analysis is 
relatively expensive in any case because it is really "piecewise 
linear;" that is, the load is applied in many small increments. After 
each load increment the new deflected shape is used to reformulate the 



Figure 2. Simply Supported, Rectangular Plate Subjected to 
a Uniform Normal Pressure Load 
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o/b =-■ 3 o 7 b = 4 

Figure 3. Finite Element Models of Rectangular Plates 


stiffness matrix. The load increments must be small enough so that 
the linear analysis within each load increment properly accounts for 
the nonlinar stiffening effect resulting from the change in geometry 
of the deflected plate. If the total load range considered is high 
relative to the load at which the behavior of the plate becomes 
nonlinear, then a large number of load increments may be required, 
driving the cost upward. A broad load range in conjunction with a 
finely divided model, which increases the cost per load increment, may 
therefore make the analysis extremely expensive. 
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The finite element models shown in Figure 3 yield reasonable 
accuracy at an affordable cost. Note that the grids shown are more 
finely divided at the corners of the plates where the stress gradient 
is highest. Even this improvement in the models incurred a cost 
penalty of the order of 50 percent. 

Figure 3 indicates the relative size of the plates analyzed. 

The absolute size of the plates is inconsequential, since the results 
are reported in terms of dimensionless parameters which describe the 
load intensity, center deflection, and stress intensity. 


B. DIMENSIONLESS PARAMETERS 


The dimensionless parameter describing the relative severity of 
the loading is called the M load intensity factor 11 and is defined as 


LIF 



( 2 ) 


The parameter used to "non-dimensional Lze M the center deflection is 
the ratio of center deflection to plate thickness, that is 


$ 


t 


Finally, the "stress intensity factor" is defined as 


SIF 


2 

ab t 
D 


(3) 


C. CENTER DEFLECTION VERSUS LOADING 

Figure 4 shows the present ARGUS results for the dimensionless 
center deflection, 6/t, as well as the results from linear theory and 
the nonlinear analytical and experimental work of other investi- 
gators. The solid lines and dots show the results of the present 
nonlinear finite element analysis. The horizontal distance between 
two consecutive solid dots represents the amount of load applied 
during a particular load increment. The vertical distance between two 
consecutive dots represents the deflection that occurred during that 
load increment. The quantity of solid dots shown is reasonably 
representative of the minimum number of load increments which may be 
taken to obtain satisfactorily convergent results over the wide load 
range considered. The reader is reminded that the abcissa of Figure 4 
is logarithmic, so the magnitude of the individual load increments may 
be continually increased as the total load increases. Smaller load 
increments are required for the rectangular plates of higher 
length/width ratios. 

The present ARGUS results (dots) for rectangular plates of 
length/width ratios of 3 and 4 are not considered satisfactory at 6/t 
ratios above 40. The reason for this problem is not known, but it is 
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suspected that cumulative errors due to the large number of load 
increments required (for length/width ratios of 3 and 4) have resulted 
in numerical problems. The recommended results in this region are 

the ?K C ’ d fo h llneS Which are extra P°lated from the results 
tor 0/ t less than 40. 

Excellent agreement of the present results with those predicted 
by linear theory is obtained for load intensity factor (LIF) values 
less than 50. In the nonlinear range, the current results for the 
center deflection of square plates are compared to the analytical 

Suearman°f R L f Vy (Ref ^ ence 9 > * the experimental results of Bowles and 
Sugarman (Reference 6), and the analytical results of Kaiser 

(Reference 10) and Tsai and Stewart (Reference 12). The current 
resuit for the center deflection of a square plate is about 25 percent 
ig er han that of Levy at a load intensity factor (LIF) of 1280. 

Levy s solution, based on a numerical solution of von Karman's 
equations m terras of trigonometric series, would appear to be 

S'* the Present problem. Beason has pointed out, however, 
at the edge conditions imposed by Levy are not identical to those 
imposed here. That is, while, the edges of Levy’s plate were free to 
ranslate laterally, they were also constrained to remain straight 
when viewed normal to the x-y plane. This accounts for the increased 
stiffness indicated by Levy’s solution. In fact, Levy stated that 
nf *1 * solution for the center deflection at a load intensity factor 
of 1280 is about 25 percent higher than his own and that "this 
difference is probably due to the fact that Kaiser allows distortions 
of the edges of the plate." The current results for the center 
deflection of square plates are in excellent agreement with the 
experimental results of Bowles and Sugarman. It is interesting to 
note that each of the data points of Bowles and Sugarman represents 
the average center deflection at fracture reached by thirty 41-inch 
square glass plates of the same thickness. The 7 points plotted are 
tor / different thicknesses and represent a total of 210 tests. The 
agreement with the analytical and experimental results of Kaiser, and 
Tsai and Stewart is also very good. Kaiser’s analysis is based on a 
finite difference solution of the von Karman equations, while Tsai and 
Stewart employed finite element analysis. Nonlinear finite element 
analysis of the center deflection of a rectangular plate of 
length/width ratio of 2, carried out by the Bechtel Corporation 
(Reference 11) using the ANSYS structural analysis program, is in 
excellent agreement with the current results. 

The linear theory is inadequate to predict accurately the 
eflection of plates at moderate to high (100 to 100,000) load 
intensity factors (Figure 4). For a square plate at load intensity 
factor of 100,000, linear theory predicts a deflection about 25 times 
that predicted by the present nonlinear finite element analysis. 


D. STRESS VERSUS LOADING 

Figures 5 and 6 show the relationship of non-dimensional ized 
stress to load (SIF vs LIF) obtained for the center and corner, 
respectiveiy, of plates of various length/width ratios. The stress 
plotted is the maximum positive principal stress. 
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Center Stress 


Excellent convergence of the present results with those 
predicted by the linear theory is seen for load intensity factor (LIF) 
values less than 100 (Figure 5). In the nonlinear load range, the 
current results are in good agreement with the analytical and 
experimental work of Kaiser, and Tsai and Stewart for square plates. 
The current results average about 5 percent higher than the former and 
13 percent higher than the latter. The present results for square 
plates are very nearly identical with the experimental results of 
Bowles and Sugarman. Levy's results for the center stress of square 
plates are also plotted for reference, even though his edge 
constraints differ from those considered here. 

The present results for the center stress of a plate having a 
length-to-width ratio of 2 are in good agreement with Bechtel's 
analytical results and with Tsai and Stewart's analytical and 
experimental results, averaging about 6 percent lower than the former 
and 2 percent higher than the latter. For plates having a length-to- 
width ratio of 4, only the analytical and experimental results of Tsai 
and Stewart are available for comparison. Their results are virtually 
identical with the current results at the low values of load intensity 
factor considered by them. As in the case of center deflection, the 
ARGUS results (dots) for the center stress for plates of length-to- 
width ratios of 3 and 4 are not considered satisfactory above load 
intensity factors of 50,000 and 30,000 respectively. The recommended 
results in this region are shown by the dot-dash lines which are 
extrapolated from the results for lower load intensities. 


2. Corner Stress 

The ARGUS structural analysis program yields the stresses at the 
center of the elements of the finite element models shown in Figure 

3. Near the center of the plates the stress changes gradually with 
respect to position in the plate. Therefore, the stress at the center 
of the element nearest to the center of the plate was assumed to be 
the same as the stress at the center of the plate. Near the corners 
of the plate, however, the stress gradient is high and the peak stress 
in the proximity of the corner of the plates cannot be assumed to be 
the 3 ame as stress at the center of the element nearest the corner of 
the plate. The peak stress in the proxmity of the corners was 
obtained by a second order La Grangian interpolation or extrapolation 
of the stresses at the center of those three elements whose centers 
lie along a path estimated to contain the point of maximum positive 
principal stress. The maximum psoitive principal stresses so obtained 
occur at or near the corner of the plate at the upper surface of the 
plate in a very localized region which is convex upward. These 
results are plotted in Figure 6. These results are in excellent 
agreement with the Bechtel analytical results for plates of 

I ength-to-width ratio of 2, being about 6 percent lower on the 
average. The current resultr are an average of 1, 12 and 8 percent 
lower than the analytical and experimental results of Tsai and Stewart 
for length-to-width ratios of 1, 2 and 4, respectively. Again, the 
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ARGUS results for length-to-width ratios of 3 and 4 are unsatisfactory 
above load intensity factors of 30,000 and 10,000, respectively. 

Better values for the corner stress may be obtained by extrapolating 
the results for lower load intensities as shown by the dot-dash lines 
in Figure 6. Comparing Figures 5 and 6, it is seen that the corner 
stress increases linearly with load to higher values of the load 
intensity factor than does the center stress. This relatively linear 
behavior of the corner stress is to be expected, since linear beahvior 
is associated with small deflections and the deflection in the close 
vicinity of the corner of the plates is comparatively small. 


E. DESIGN CURVES FOR STRESS VERSUS LOADING 

It can be seen (Figures 3 and 6), that stress at the center of 
the plates is higher than that at the corners of the plate for low to 
moderate loading. At higher loads, however, the stress at the corners 
of the plate is higher. Figure 7 is constructed by superimposing 
Figure 5 on Figure 6, but showing only the higher of the center or 
corner stresses at any given load intensity. Figure 7 yields the 
maximum positive principal stress, a, versus load intensity factor for 
length- to-width ratios of 1, 1.5, 2, 3 and 4. The maximum positive 
principal stress occurs on the bottom surface of the plate at the 
center and on the top surface of the plate at the corners. At the 
center it is predominantly membrane tension stress plus bending 
stress. At the corners it is largely shear stress. 
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Figure 7. SIF versus LIF Showing Composite Curves of the Larger of 
the Maximum Positive Principal Stresses on Plate (Center 
Bottom or Top Surface Near Corner) 






SECTION VI 


ALLOWABLE STRESS FOR GLASS 


A brief discussion of the factors influencing the failure of 
glass, adapted from Brown (Reference 14), is included here. In 1921 
Griffith noted the weakening effect of preexisting flaws or scratches 
in the glass surface. The high stresses at the tips of these flaws 
account for the fact that the measured strength of glass in tension is 
usually several orders of magnitude lower than its theoretical 
strength. Deliberate surface conditioning of the glass surface by 
various investigators has borne out Griffith's observation. Surface 
roughening reduces the mean strength and variability. Polishing 
increases the mean strength and variability. The measured strength of 
glass samples tested in moist air decreases with increased duration of 
loading. It is now generally accepted that slow flaw growth takes 
place as a result of stress corrosion. The mathematical model 
describing the dependence of glass failure on load duration and 
loading proposed by Brown (Reference 14) has been adopted by Beason 
(Reference 4) and Dalgliesh (Reference 3) in the following form. 



= constant 


(4) 


where 


T ■ t ime 

Tg * time to break 

a(l) * applied, time-variant stress 

a * an empirically determined constant which depends on the 
surface condition of the glass, relative humidity and 
temperature . 

With these factors considered, a quantitative assessment of the 
effect of glars temper, time duration of loading, glass surface area, 
and probability of failure is presented in the following paragraphs. 


A. EFFECT OF TEMPER AND LOAD DURATION 

As discussed above, the breakage strength of glass decreases as 
the duration of loading increases. For a constant applied stress 
Equation 4 may be written in the form 


1 
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where 

Og * breakage strength 

t b = time to break 

a 3 an empirically determined constant, applicable only over a 
narrow range of load durations where it fits the data. 

Dalgliesh (Reference 15) reports values of a for annealed glass 
from 12 to 20. Dalgliesh does not recommend a value of a for tempered 
glass. Shand (Reference 16) states that not only is tempered glass 
stronger initially, but the decrease in the breaking strength of 
tempered glass with load time is significantly less than that for 
annealed glass. This implies that the value of a would be signifi- 
cantly higher than 20. 

For the simple power law (Equation 5) just discussed, the 
fracture 9tress approaches zero for loads of very long duration. This 
is at odds with the widely accepted idea that there is a level of 
applied stress, called the "endurance limit, 11 below which crack growth 
will not occur. At levels of applied stress less than the endurance 
limit glass parts should be able to sustain loads for an indefinite 
period of time. Figure 3, reproduced from Shand (Reference 16), shows 
the breaking stress versus the duration of that stress. These curves 
are consistent with the concept of an endurance limit in that they 
tend toward some asymptotic lower bound of the breaking stress for 
long-term loading. Unfortunately, these curves extend only to a load 
duration of 10? seconds (approximately 4 months), whereas to 
properly consider deadweight load, we need load durations up to the 
design life of the glass plates (that is, greater than 20 years). 

Also, it should be noted that Shand (Reference 16) does not specify 
the size of the glass specimens on which Figure 8 is based. 



Figure 8. Stress-Time Characteristics of Class Broken in Flexure 
Tests at Room Temperature (Composite Curves) 
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Therefore, it is the relative strength of annealed glass versus 
tempered glass and the decrease in strength with time that are of 
interest here. 

The curves shown in Figure 9 have been extrapolated to a load 
duration of approximately 100 years and normalized to a value of unity 
for the 1-minute breaking stress of the glass. Shand (Reference 16) 
has stated that the endurance limit for (annealed*) glass is of the 
order of 40 to 45 percent of the 5-second breaking strength. The 
curve for annealed glass (Figure 9) has been faired to a value of 
about 40 percent of the 5-second breaking strength at 100 years load 
duration. This is somewhat less than if Shand's original curve (see 
dotted line - Figure 9) had been simply extrapolated to longer load 
durations. Likewise the curve for tempered glass shown in Figure 9 is 
somewhat lower at long load durations than Shand* s original curve 
(shown dotted). To obtain strength-versus-t irae values for 
semi-tempered glass, the values \or annealed and tempered glass have, 
arbitrarily, been averaged. The resulting values for semi-tempered 
glass are shown as a dot-dash curve in Figure 9. The equation 
f T 7 ( V 1 ) 1 ' 1 * 2 is also plotted in Figure 9. This equation is 
equivalent to Equation 5 and shows that for load durations between 5 
seconds and 1 minute the value of a equal to 12 provides an excellent 
fit to the annealed glass curve. Likewise, a value for a of 40 
provides an excellent fit to the tempered glass curve for load 
durations between 5 seconds and 5 minutes. 


B. EFFECT OF PLATE AREA 

As mentioned earlier, large glass parts fail at lower applied 
stress levels than do small glass parts. This is true because the 
likelihood that a flaw in the glass will coincide with a region of 
high applied stress is greater for larger parts. These flaws exist in 
newly manufactured glass and their size and/or number increase with 
time due to handling, missile impact and other loads. In the vast 
majority of cases, failure originates from surface flaws so that the 
area of a glass plate represents the pertinent measure of part. size. 
Dalgliesh (Reference 15) reported that laboratory tests on glass 
plates have shown that the breaking strength varies inversely as the 
fifth to seventh root of glass surface area. 


1 



*It is assumed that Shand means annealed glass since he further states 
that the decrease in the breaking strength of tempered glass with 
increasing load duration time is significantly less than that for 
annealed glass. 
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Figure 9. Fraction of 1-Minute Load Duration Glass Breakage 
Stress versus Load Duration 


For the purpose of the design method presented herein, it is 
expedient to define the fraction, f A , which is the fraction of the 
breaking stress of a 1-square-meter plate which will be attained by a 
plate of area A if both plates break 1 minute after the sudden 
application of the full load. From Equation 5 the fraction f A may 
be expressed mathematically as 
1 


where 

f. = fraction of the breaking stress of a 1 square meter 
plate which will be attained by a plate of area A 

A = area of plate for which breaking stress is unknown 
(expressed in square meters) 


Equation 7 is plotted on Figure 10. 
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Figure 10, Fraction of Breakage Strength of 1-Squaro 
Motor Class Pintos versus Plato Aroa 


0. BREAKAGE STRENGTH VERSUS PROBABILITY OF FAILURE 

In the following paragraphs the glass plate breakage data of 
other investigators is reanalyzed to obtain the glass breakage 
strength as a function of the probability of failure. 


1 . Selocti on of G 1 as s P 1 a t: e Brea kage Da t a 

To obtain a practical value for the breaking strength of glass 
to he used in determining the required thickness of photovoltaic sola 
panels the raw data of Bowles and Sugarman (Reference 6), Beason 
(Reference B), and Wilson (Reference 7) have been analyzed employing 
Wei. bull statistical analysis techniques. The data of these 
investigators are preferred for the following reasons: 

( a ) The size o f the s am pies tested is i n t e rme d i a t e be t we en t b 
smallest and largest glass photovoltaic modules being 
considered. Bowles and Sugarman tested 41-inch square 
glass plates, Beason tested 28.^ x 60. S- inch plates and 
28,‘j-inch square plates. Wilson tested 48-inch square 







(b) These investigators employed edge constraints which 
approximate a simply-supported condition. 

(c) The load-time history prior to panel failure is necessary 
to make a meaningful assessment of the breaking stress for 
other load-time histories. The load-time histories from 
these data sources are known reasonably well. Bowles and 
Sugarman loaded the panels with a pressure load which 
increased approximately as the second power of time. They 
adjusted the loading mechanism so that the average 
pressure resulting in failure was reached 30 seconds after 
the start of the test. Reason 1 s raw data are more 
definitive in this respect. His data include a complete 
pressure-time history for each sample of glass tested, as 
does Wilson's. 

(d) These investigators tested a sufficient number of samples 
for the results to be statistically significant. Bowles 
and Sugarman tested 40 samples of .122-inch-thick annealed 
plate glass, 30 samples each of .197, .250, .373-inch- 
thick annealed plate glass and .110, .158, .195-inch-thick 
annealed sheet glass - all 41 inches square. Beason 
tested 20 samples each of .219 x 28.5 x 60.5-inch and .219 
x 28.5 x 28.5-inch annealed sheet glass. In Beason's 
tests the 20 samples of each size were divided into 10 
samples with the weathered side in tension and 10 samples 
with the indoor side in tension. Wilson tested 8 samples 
of .125 x 48 x 48-inch tempered float glass. 

(e) The Bowles and Sugarman tests were performed on new 
annealed sheet and plate glass as noted above. Beason's 
tests were performed on "weathered" glass, removed from 
the Great Plains Life Building, Lubbock, Texas. This 
glass had been in service for about 20 years at the time 
of its removal for testing; it is believed to be annealed 
sheet glass. Wilson tested new tempered float glass. 

Tables 1 and 2 give the burst pressure (p^) data for the 
Bowles and- Sugarman tests of annealed plate glass and annealed sheet 
glass, respectively. These data were adapted from his tograramic data 
included in the Bowles and Sugarman report. Tables 3, 4, 5 and 6 show 
the pressure-time history of Beason's tests of the 20-year-old 
weathered annealed sheet glass plates. Table 7 shows the 
pressure-time history of Wilson's tests of new tempered float glass. 


2. Normalizing the Glass Plate Breakage Data 

Examination of the data in Tables 1 through 7 shows that glass 
plates of various sizes break over a wide range of pressures after 
experiencing different pressure-time histories. It has already been 
noted that the breakage strength of glass is a function of the load 
duration and the surface area of the glass plate. 
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Table 2. Data for New Annealed Sheet Glass 


Burst pressure data for 41-inch square, simply supported, new 
annealed sheet glass plates (from Bowles and Sugarman - 
Reference 6) 


t 

Thickness 

(in) 

Number of 
Panels 
Failed 

P B Burst 
Pressure 
(lb/ in^) 

t 

Thickness 

(in) 

Number of 
Panels 
Failed 

P B Burst 
Pressure 
(lb/ in^) 

0.110 

1 

0.35 

0.195 

i 

1.05 


i 

0.55 


i 

1.15 


13 

0.65 


i 

1.25 


11 

0.75 


i 

1.35 


4 

0.85 


1 

1.55 





3 

1.65 

0.158 

3 

0.95 


3 

1.75 


2 

1.15 


2 

1.85 


5 

1.25 


3 

1.95 


8 

1.35 


4 

2.05 


4 

1.45 


2 

2.15 


2 

1.55 


2 

2.25 


5 

1.65 


2 

2.35 


1 

1.75 


3 

2.45 





1 

2.55 



Table 3 . Data for Weathered Annealed Sheet Glass 
(weathered side in tension) 


Pressure-time history data for 0.219 x 28.5 x JO.S-ineh simpiy 
suoDorted, 20-year-old weathered annealed sheet glass plates 
fiom Great Plains Uf. Building, bubbock, Texas (fro. 

Beason - Reference 5 ). 


Sample 

No. 


T 

P 

Time 

Pressure 

(sec) 

(lb/ in 2 ) 

5.0 

0.090 

10.0 

0.20b 

15.0 

0.344 

20.0 

0.449 

25.0 

0.548 

30.0 

0.626 

33.5 

0.664 

5.0 

0.141 

10.0 

0.256 

15.0 

0.364 

20.0 

0.475 

25.0 

0.568 

30.0 

0.641 

35.0 

0.710 

40.0 

0.760 

42.0 

0.765 

5.0 

0.112 

10,0 

0.225 

15.0 

0.335 

20.0 

0.430 

25.0 

0.515 

30.0 

0.602 

35.0 

0.691 

40.0 

0.771 

45.0 

0.867 

50.0 

0.939 

55.0 

0.986 

56.7 

0.999 

5.0 

0.090 

10.0 

0.191 

15.0 

0.304 

20.0 

0.420 

25.0 

0.520 

30.0 

0.600 

35.0 

0.686 

40.0 

0.737 

41.5 

0.742 

5.0 

0.093 

10.0 

0.179 

15.0 

0.269 

20.0 

0.377 

25.0 

0.456 

30.0 

0.518 

35.0 

0.579 

40.0 

0.609 

45.0 

0.653 

50.0 

0.668 


Sample 

No. 


Pressure 
( lb/in 2 ) 
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Table 4 


Data for Weathered Annealed Sheet Glass 
(weathered side in tension) 


Pressure-time history data for 0.219 x 28.5 x 28.5-inch simply 
supported, 20-year-old weathered annealed sheet glass plates 
removed from Great Plains Life Building, Lubbock, Texas (from 
Beason - Reference 5). 
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Table 5 


. Data for Weathered Annealed Sheet Glass 
(indoor side in tension) 


Pressure-time history data for 0.219 x 28.5 x 60.5— inch simply 
supported, 20-year-old weathered annealed sheet glass plates 
removed from Great Plains Life Building, Lubbock, Texas (from 
Beason - Reference 5). 


Sample 

No. 

Time 
( sec ) 

P 

Pressure 
( lb/ in^) 

Sample 

No. 

i 

T ime 
(sec ) 

P 

Pressure 
( lb/in^) 

3 

3.0 

0 , !1 004 

15 

5.0 

0.095 


8.0 

0.122 


10.0 

0.217 


13.0 

0.269 


15.0 

0.318 


18.0 

0.379 


20.0 

0.416 


23.0 

0.446 


21.2 

0.449 


28.2 

0.535 







17 

5.0 

0.104 

7 

3.0 

0.077 


10.0 

0.273 


10.0 

0.168 


15.0 

0.434 


13.0 

0.267 


18.0 

0.530 


20.0 

0.372 





24.7 

0.462 

19 

5.0 

0.086 





10.0 

0.144 

9 

5.0 

0.097 


15.0 

0.250 


10.0 

0.207 


20.0 

0.339 


15.0 

0.29.5 


25.0 

0.425 


20.0 

0.368 


30.0 

0.502 


25.0 

0.451 


35.0 

0.550 


30.0 

0.542 


40.0 

0.587 


35.0 

0.630 


44.7 

0.611 


40.0 

0.721 





45.0 

0.807 

21 

5.0 

0.104 


50.0 

0.888 


10.0 

0.225 


51.2 

0.901 


15.0 

0.374 





20.0 

0.503 

11 

5.0 

0.069 


25.0 

0.604 


10.0 

0.153 


28.8 

0.670 


13.0 

0.222 





20.0 

0.267 

23 

5.0 

0.102 


23.0 

0.290 


10.0 

0.221 


25.0 

0.346 


15.0 

0.314 


30.0 

0.442 


20.0 

0.404 


35.0 

0.514 


25.0 

0.457 


36.2 

0.532 


30.0 

0.492 





34.5 

0.516 

13 

5.0 

0.127 





10.0 

0.255 





15.0 

0.350 





20.0 

0.423 





25.0 

0.482 





30.0 

0.516 





35.0 

0.547 





40.0 

0.556 





42.5 

0.562 









immmuw 


Table 6. Data for Weathered Annealed Sheet Glass 
(indoor side in tension) 

Pressure-time history data for 0.219 x 28.5 x 28.5-inch simply 
supported, 20-year-old weathered annealed sheet glass plates 
removed from Great Plains Life Building, Lubbock, Texas (from 
Beason - Reference 5). 



T 

P 


X 

P 

Sample 

Time 

Pressure 

Sample 

Time 

Pressure 

No. 

(sec) 

(lb/in^) 

No • 

(sec) 

(lb/ in 2 ) 

18 

9.50 

o.^o 

23 

3.30 

0.137 


19.50 

0.716 


13.46 

0.455 


29.50 

1.046 


23.62 

0.849 


39.50 

1.401 


27.08 

0.955 

19 

4.25 

0.156 

24 

5.08 

0.273 


14.25 

0.484 


15.28 

0.636 


24.25 

0.835 


25.40 

0.985 


34.25 

0.124 


35.56 

1.288 





37.59 

1.363 

20 

6.5 

0.258 





16.5 

0.571 

25 

10.16 

0.332 


26.5 

0.939 


20.32 

0.740 


36.5 

1.229 


30.48 

1.149 


46.5 

1.621 


35.05 

1.300 

21 

4.06 

0.166 

26 

3.64 

0.166 


14.22 

0.529 


13.80 

0.513 


24.38 

0.878 


23.96 

0.906 


34.54 

1.197 


25.59 

0.997 


42.21 

1.438 

27 

1.78 

0.135 

22 

4.32 

0.196 


11.94 

0.453 


14.48 

0.604 


22.10 

0.816 


24.64 

0.920 


32.26 

1.073 


34.08 

1.176 


42.42 

1.632 


37.74 

1.313 


43.94 

1.647 
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Table 7 


Data for New Tempered Float Glass 


Pressure-time history data for 0.125 x 48 x 48-inch simply 
supported, new tempered float glass plates (from Wilson — 
Reference 7). 
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The data in Tables 1 through 7 are normalized to a breaking stress 
corresponding to a constant stress applied to a plate of 1 square meter 
surface area for 1 minute. This normalization is accomplished as 
follows: 

(a) First, the pressure-time histories for the various tests 
are converted to stress-time histories, using Figure 7 and 
Equations 2 and 3. This is necessary since the stress-time 
history is needed and the stress-pressure relationship is 
nonlinear. 

8 

(b) An analytical expression of the form a 9 ki is fitted to 
the above stress-time data. 


(c) 


This stress-time relationship together with Equation 4 
(repeated here) 



* constant 


(4) 


is used to determine aj, that constant stress which vill 
cause the panel to break after 1-minute load duration. 

For a glass load-bearing plate which fails at time Tjj 
and stress cr B , to which the stress is applied gradually 
according to the relationship a = ki^, and bearing in mind 
that we want to develop a relationship to normalize the 
data to the 1-minute breaking stress, (J i , where i is in 
seconds and Tg = 60 seconds, we substitute in Equation 4 


f 60 a f l i 

Jo “'’Ho 


8 o nt 
( k v ) dT 


Simplifying, we obtain 


° l ' [60(,«6. l)] 1/a 
where 


( 8 ) 


(j j » constant stress that will cause the panel to 

break at 1-rainute load duration 

k,S • empirical constants obtained by fitting an 

analytical expression to the stress-time data of 
step b above. 

a ■ empirical constant associated with the decrease 
of glass breakage stress with increased loaa 
duration. From Figure 9 we note that values of a 
equal to 12 for annealed glass and to 40 for 
tempered glass fit the data for load durations in 
the 5-second to 5-minute range. 
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Tjj = time in seconds at which a given glass panel 
fails. 


(d) Note that the 1-minute breaking stress obtained for each 
plate tested applies only to plates having the same 
surface area as the test plate. Employing Equation 6, the 
one minute breaking strengths are further normalized to 
the 1-minute breaking strength corresponding to a glass 
plate with a surface area of 1 square meter, as follows: 


a 


11 


1/6 

= (A) a x 


(9) 


where 

Ojj * constantly applied stress that will cause a panel 
of 1- square-meter surface area to break at 
1-minute load duration. 

A * surface area of the test plate expressed in 
square meters. 

o, a constantly applied stress that will cause the 
test panel of A square meters to break at 1 
minute load duration. 

The normalized breaking strengths given in Tables 8 through 14 are 
produced following the procedure outlined above. 


3. Breaking Strength versus Probability of Failure 

Weibull statistical analysis techniques are applied to the 
breakage stress data of Tables 8 through 14 to obtain analytical 
expressions for the probability of failure as a function of the 
breakage stress. The use of the Weibull distribution represents 
state-of-the-art statistical analysis of the fracture of brittle 
materials. The Weibull distribution generally provides a better fit 
to experimental fracture data than does the Gaussian distribution, 
especially at low probabilities of failure. 
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Table 8. Normalized Breakage Stress for New Annealed Plate Glass 


Bowles and Sugarman's burst pressure data for 41-inch square, simply 
supported, new annealed plate glass plates (Table 1) normalized to 
°11» that constantly applied stress that will cause a plate of 
1-square-meter surface area to fail at 1-minute load duration. 



Number 

°11 

1-Minute 


Number 

°11 

1-Minute 

t 

of 

Load Duration 

t 

of 

Load Duration 

Thickness 

Plates 

Stress for 1 

Thickness 

Plates 

Stress for 1 m 2 

(in) 

Failed 

Plate (lb/in 2 ) 

(in) 

Failed 

Plate (lb/in 2 ) K 

0.122 

i 

4083 

0.373 

1 

4569 


5 

4802 


1 

4849 


3 

5496 


1 

4984 


13 

6174 


1 

5120 


12 

6834 


2 

5253 


4 

7475 


1 

5515 


1 

8109 


2 

5646 


1 

8728 


2 

5902 





3 

6027 

0.197 

1 

4301 


1 

6275 


2 

4675 


2 

6396 


5 

5039 


2 

6517 


3 

5401 


3 

6757 


3 

5754 


1 

6877 


4 

6106 


I 

7111 


3 

6452 


2 

7458 


4 

6795 


2 

7684 


4 

7135 


1 

8454 


1 

7470 


1 

8776 

0.250 

1 

3683 





3 

3956 





4 

4225 





1 

4490 





1 

4754 





2 

5013 





3 

5278 





7 

5529 





2 

5784 





1 

6536 





1 

6781 





1 

7027 





1 

7274 





1 

7757 





1 

7995 
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Table 9. Normalized Breakage Stress for New Annealed Sheet Glass 


Bowles and Sugarman's burst pressure data for 41-inch square, simply 
supported, new annealed sheet glass plates (Table 2) normalized to 
o n , that constantly applied stress that will cause a plate of 
1-square-meter surface area to fail at 1-rainute load duration. 


Number 

t of 

°11 

1-Minute 
Load Duration 

t 

Number 

of 

°n 

l-Minube 
Load Duration 

Thickness 

Plates 

Stress for 1 m^ 

Thickness 

Plates 

Stress for 1 m 
Plate (lb/in^) 

fin) 

Failed 

Plate (lb/in^) 

(in) 

Failed 






Table 10- Normalized Breakage Stress for 
Weathered Annealed Sheet Glass 
(weathered side in tension) 

Season’s burst pressure data for 0.21.9 x 28.5 x 60.5-inch, simply 
supported, 20-year-old weathered annealed sheet glass removed from 
Great Plains Life Building, Lubbock, Texas (Table 3) normalized to 
° 1 1 , that constantly applied stress that will cause a plate of 
1-square-meter surface are.* to fail at 1-minute load duration. 


t 

Thickness 

(in) 


Number of 
Plates Failed 


1-Minute Load Duration Stress 
for 1 m^ Plate (lb/in^) 


0.219 1 4438 

1 4562 

1 4934 

1 5052 

1 5079 

1 5285 

1 5370 

1 5421 

1 5700 

1 6524 


Table 11. Normalized Breakage Stress for 
Weathered Annealed Sheet Glass 
(weathered side in tension) 

Season’s burst pressure data for 0.219 x 28.5 x 28.5-inch, simply 
supported, 20-year-old weathered annealed sheet glass removed from 
Great Plains Life Building, Lubbock, Texas (Table 4) normalized to 
^\\, that constantly applied stress that wilt cause a plate of 
1-square-meter surface area to fail at l -minute load duration. 


t 

Thickness 

(in) 


Number of 
Plates Failed 


1 -Minute Load Duration Stress 
for l m*~ Plate (lb/in^) 


0.219 l 2309 

l 2690 

l 3306 

1 3471 

l 3801 

1 3894 

1 4094 

t 4127 

l 4573 

l 5188 
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i 


Table I.’* Normalized Breakage Stress tor 
Weathered Annealed Shoot Class 
(indoor si do in tons ion) 


Be a son 1 s hurst pressure data tor 0.219 x 28. > x bO.h-iiuht simply 
supported, 20-vear-old woathorod anuoalod shoot glass removed from 
Croat Plains l.it'o Building, l.ubboek, Texas (Table r >) normalized 
that constantly applied stress that will cause a plate of 
l -square-mot or surface area to tail at l -minute load duration. 


to 


t 

Th i ok ness 

( in ) 


Number of 
Plates Failed 


l ’ll 

l-Minuto Load Duration Stress 
f o r l P late ( l h/ i n - ) 


O.'MO 1 

1 

1 

l 

1 

1 

l 

1 

1 

l 


1200 

3312 

3S79 
3802 
38 tS 
3 i )2(> 
A3(>7 
AS10 
h S 2 3 
S804 


Table 13. Normalised breakage Stress for 
Weathered Annealed Sheet Class 
(Indoor Si vie in Tension) 


Reason's burst pressure data tor 0.719 x 78.^ x .8.'- 
avipported, 20-vea r-o Id weathered annealed sheet glass 
Croat Plains Lite lluilvlin.o., l.uhhook, Toxas (TaMo 6) 
•'j), that oonstantlv appliovl st ross that will ‘aviso a 
l-squaro-mot or surfaoo aroa to fail at l-minuto load 


- i noli , simply 
removed from 
normalised to 
p 1 a t o o f 
dura t ion. 


t 

Th i cknes s 
( in) 


Number of 
Plates Failed 


v \ l 

1-Minute l.oad Duration Stress 

for l m* Plate (lb/ in*) 


0.210 


l 

l 

l 

i 

i 

i 

t 

i 

i 

i 


2662 

2741 

1088 

14 It 

15 07 
IhOl 
l6l‘> 
1781 
1087 
407 2 


6 - 1 ') 


M 


* 


Table 14. Normalized Breakage Stress for 
New Tempered Float Glass 

Wilson's burst pressure data for 0.125 x 48 x 48-inch, simply 
supported, new tempered float glass (Table 7) normalized to 

that constancly applied stress that will cause a plate of 
1-square-meter surface area to fail at 1-minute load duration. 


t 

Thickness 

(in) 

Number of 
Plates Failed 

CT n 

1-Minute Load Duration Stress 
for 1 Plate (lb/in^) 

0.125 

i 

16986 


i 

17344 


i 

17816 


i 

17966 


i 

18031 


i 

18531 


i 

19264 


i 

19766 


The analytical expression for the Weibull distribution used here is 



where 

P f = probability of failure 

= breakage stress for 1-minute load duration for 
1-square-meter plate 

a u * lower bound stress for which the probability of failure 
is zero. 

C Q * characteristic stress 
m ■ Weibull modulus 


The constants ci u , o 0 and m for a given data set are determined 
from a curve-fitting method described by Sines (Reference 
basic procedure is as follows; 

( a ) Take the logarithm of both sides of Equation 10 twice to 
produce the following relationship: 

log log ( r ~jT> = ® log <o U " lQ g " 0.362216 

(b) For any given data set of Tables 8 through 1*+ assume 
value of 

(c) compute the experimental P f associated with each value 
of & \ i : 

>’f ■ FTT 


(d) 

(e) 

(f) 

(fi> 

(h) 


where 

n = cumulative number of failures at or below a given < T n . 

N = total number of samples in a given data set. 

Now plot log log ( T 4 t 7 ) VCrSUS l ° g ( °11 " CT u ) imd CO “ PUte 
the best fit straight line through these points by the 
method of "least squares." 

Compute the correlation coefficient which measures '■»” uel > 
the above heat fit straight line fits the experimental 
data. The correlation coefficient will have a value between 
zero and one, with one indicating a perfect fit. 

Varv it,, until the correlation coefficient attains a 
maximum value (that is, closest to one), tht, is he va.uo 
of o„ for the Heibull distribution which best fits the 
experimental data* 

The Wcibull modulus, m, is the slope of the best fit 
straight line corresponding to the above value of o u . 

Having a u and m, compute a 0 from the relationship of 
step (a). We now have all of the constants (a a n , m) 
to form the Weibull distribution (Equation 10) which bos 
describes a given data set. 
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Robert Weaver of Jet Propulsion Laboratory wrote a FORTRAN 
computer program to automate the above procedure. A listing of this 
program together with a sample output is given in Appendix II. This 
analysis was performed on each of the 12 normalized data sets of Tables 
8 through 14. Figure 11 shows the resulting "best fit 11 Weibull 
distributions (lines) together with the experimental Pf, pairs 

(symbols) . 

In comparing the analytical Weibull distributions for the various 
annealed glass types (new plate, new sheet and weathered sheet), it is 
very interesting to note that the mean strength of the new annealed 
sheet glass samples is generally higher than that of the new annealed 
plate glass samples. Moreover, the lower bound strength (p£ » 0) of 
the new annealed plate glass samples has some appreciable value (a u 
greater than about 3200 lb/in^) while the lower bound strength of all 
new annealed sheet glass samples is zero. This can be explained by 
considering the probable nature of the flaws in the surface of these 
glass types resulting from the different manufacturing methods used to 
produce plate and sheet glass. 

Sheet glass is produced by drawing the molten glass from a bath. 
The process is continuous and the glass cools as it is drawn over 
rollers. No subsequent operations are performed. The general surface 
condition of sheet glass might be described as " f ire-pol ished" , which 
results in the high mean strength noted for sheet glass. On the other 
hand, the chance for a large surface inclusion or flaw exists, which 
explains why the new annealed sheet glass samples have zero lower bound 
strength. The final operation in the manufacture of plate glass is 
grinding the surface smooth and flat. This results in a fairly even 
surface flaw size distribution, with the large inclusions and flaws 
having been removed. This accounts for the lower mean strength and the 
substantial lower bound strength of plate glass. The strong dependence 
of the measured mean strength and strength distribution on the surface 
condition of the glass has been noted by others. Brown (Reference 14) 
states: "Surface roughening results in reduced average strength and 

reduced variability, whereas smoothing operations result in increased 
strength with attendant increased variability." 

It is interesting to note the different character of the Weibull 
d i s tr ibut ions * for the four lots of weathered annealed glass tested by 
Beason. Beason tested two lots of 12 plates each at the original size 
removed from the Great Plains Life Building ir. Lubbock, Texas; that is, 
0.219 x 28.5 x 60.5 inches. The additional two lots tested included 10 
samples each of plates 0.219 x 28.5 x 28.5 inches cut from the original 
plates. For both sizes, 10 plates were tested with the weathered side 
in tension and 10 plates with the indoor side in tension. The data from 
the plates tested with the indoor side in tension are included for 
comparative purposes and because they are generally supportive of the 
following conclusions: 


LASS BREAKAGE STRENGTH ~ irsi 
-m 2 PLATES. 1 -min LOAD DUPATi 





(a) Comparing the Weibull curves for the two different plate 
sizes tested with the indoor side in tension supports the 
use of the 1/6 power areal correction factor, 

(b) Since both sample lots of 28.5 x 28.5-inch plates have zero 
lower bound strength (a u - o) whereas both sample lots of 
28.5 x 60.5-inch plates have substantial lower bound 
strength, cutting the glass presumably resulted in 
excessive surface or edge damage. Therefore, the data from 
the 0.219 x 28.5 x 28.5-inch samples is invalid, at least 
for assessing the strength of naturally weathered glass. 

(c) The only sample lot which is valid for assessing the 
strength of weathered annealed glass is the one containing 
ten 0.219 x 28.5 x 60.5-inch plates tested with the 
weathered side in tension. 

Figure 12 shows recommended design values for breaking strength 
versus probability of failure for new annealed plate glass, new 
annealed sheet and float glass, and weathered annealed glass. It was 
derived from Figure 11 as follows. For new annealed plate glass the 
mean of the Figure 11 curves for thicknesses of 0.122 inch, 0.197 inch 
and 0.373 inch is plotted in Figure 12. The curve for 0,250-inch 
thickness is not included in this mean as it is sufficiently different 
from the other curves to be nonrepresentative of new annealed plate 
glass. For new annealed sheet glass the mean of all three thicknesses 
is plotted in Figure 12. As discussed previously, only the curve shown 
in Figure 11 for the 0.219 x 28.5 x 60.5-inch weathered annealed sheet 
glass samples tested with weathered side in tension is taken as 
representative of weathered annealed sheet glass and is reproduced in 
Figure 12. The curve shown in Figure 11 for 0.125 x 48 x 48-inch 
tempered float glass is reproduced in Figure 12. 

The curve for annealed sheet glass is recommended for annealed 
float glass also, as the processing would most likely produce the same 
type of surface flaws. The curve for weathered annealed sheet glass is 
recommended for all types of weathered annealed glass, including plate, 
sheet and float glass. At a 1Z* probability of failure the total range 
of breaking strengths for annealed glass plates, regardless of glass 
type or age, varies only from 3800 lb/in^ to 4000 lb/in^. 

The curve for 0.125-inch-thick tempered glass should not be used 
for thinner tempered glass because of the difficulty of achieving a 
high level of initial surface compressive stress in glass plates 
thinner than 0.125 inch. On the other hand the strength of tempered 
glass plates thicker than 0.125 inch would be somewhat higher than this 
curve . 


♦Architects generally employ a design breakage rate of 8 per 1000, 

(P f - 0.8 Z). 


LASS BREAKAGE STRENGTH - ksi 
•rr, 2 Plates, 1-mln LOAD DURATION 



Figure 12. Recommended Design Values for Breakage Strength 
versus Probability of Failure for 1 m 2 , Simply 
Supported, Glass Plates Subjected to a Uniform 
Normal Pressure Load of 1-Minute Duration 



SECTION VII 


METHODOLOGY 


A rational means of determining the required thickness of a 
simply supported, rectangular glass plate subjected to environmental 
loads which may be characterized as uniform normal pressure loads is 
described in detail in the following pages* The method is outlined 
below: 

(1) Define the various loading conditions (wind, earthquake, 
snow, deadweight, etc.) appropriate for a given locale 
using local building codes, or the codes of other advisory 
or regulatory agencies, such as ANSI A58.1 - 1972 
(Reference 18), as applicable. 

(2) Determine the stress level corresponding to each of the 
above loading conditions for an assumed glass thickness. 

(3) Determine the allowable stress levels corresponding to 
each of the above loading conditions. The allowable 
stresses are a function of the required failure rate, 
degree of tempering of the glass, surface area of the 
glass, and dur*‘ .1 of the load. 

(4) Compare the estimated stresses to the allowable stresses 
for various load combinations to see if the assumed glass 
thickness is adequate. 


A. DEFINE THE LCADINC 


Define the uniform normal pressures associated with wind, snow, 
deadweight and other loads. The following procedure should be 
employed : 


(l) Using ANSI A58.1 - 1972 (Reference 18) or another code 
appropriate for the locale, define the various loading 
conditions in terms of a uniform normal pressure applied 
to the plate. Thermal loads should also be included if 
they produce a tensile stress in the glass plate. The 
designer should consider at least the following loads. 


Load Description 


Pressure 


Wind 

Earthquake 

Snow 

Deadweight 

Thermal 


(Tensile stress in 
glass obtained from 
other analysis) 
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(2) If the glass plate being analyzed is the only structural 
member which resists the above loads, skip this section. 
If, however, the glass plate being analyzed is part of a 
sandwich structure, it is now necessary to calculate that 
portion of the load which the glass plate will bear. It 
should be noted that the normal encapsulant/adhesives 
typically employed in the manufacture of photovoltaic 
solar panels do not have sufficient sheer stiffness to 
cause the top and bottom faces of these sandwich struc- 
tures to act together as a composite plate. Therefore, 
the only structural benefit to be derived from such sand- 
wich configurations is that the pressure loads are shared 
according to the relative stiffness (flexural regidity, D) 
of the glass plate being analyzed and the other laver(s) 
of the sandwich. Therefore, the pressure to be applied to 
the glass plate may be calculated as follows. A typical 
sandwich-type photovoltaic solar panel is shown in Figure 
13 to help illustrate the method of apportioning the 
pressure loads among the layers of the sandwich. 



Figure 13. Cross-Section through Typical 

Sandwich-Type Photovoltaic Panel 


Calculate the flexural rigidities of the various layers: 


Layer 


Flexural Rigidity 


Glass 


Encapsulant/ 

adhesive 


Additional load- 
sharing layers 



12(l-v ) 



12(1 - v' ) 


D M 


E'V' 3 

12( 1- v" 2 ) 


(The stiffness 
of this layer 
may normally 
be neglected) 


Determine the appropriate fraction of the pressure loads 
of step Al, p'/pi that should be applied to the glass 
layer being analyzed. 

El D 

p D + D' + D" + . . . 
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1 


Modify the table of loads from step Al to obtain the 
pressures borne by the g^3ss layer being analy eds 


Load Description 


Pressure 


Wind 

Earthquake 

Snow 

Deadweight 

Thermal 



(Tensile stress in 
glass obtained from 
other analysis) 


B. ESTIMATE THE STRESSES 


Now determine the stress levels corresponding to each of the 
above loading conditions. 

(1) Assume a glass thickness, t. 

(2) Calculate the "Load Intensity Factor" (LIF) corresponding 
to each of the above loading conditions. The 
dimensionless LIF is a function of the loading, and the 
geometry and material properties of the plate as follows: 

LIF » (2) 

where 

p = uniform normal pressure 
b M width of plate 
t = thickness of plate 
D = flexural rigidity of plate 



12(l-v 2 ) 


E * Young's modulus of plate (use 10,000,000 lb/in^ 
for glass) 

- Poisson's ratio of plate (use 0.22 for glass) 

(3) From Figure 7 determine the "Stress Intensity Factor" 

(SIF) corresponding to the Load Intensity Factor (LIF) for 
each of the loading conditions of step Al. Note that the 
SIF is a dimensionless quantity related to the maximum 
positive principal stress which will occur in the plate. 
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a 


The maximum value may occur in either the center of the plate 

or the corner of the plate depending on the intensity of the 
loading. 


(4) The maximum positive principal stress corresponding to each 

of the loading conditions is now calculated from the equation 


(5) 


a 



(3) 


Now construct a table of stresses corresponding to the 
various loading conditions 


Loading 

St 

Wind 


Earthquake 

°\ 

Snow 


Deadweight 


Thermal* 

i 

a, 


C. 


DETERMINE THE ALLOWABLE STRESS FOR GLASS 


The allowable stress for a given loading condition is the 

cIr^?!/ t f r ' S ^ tr ° m ! igure 12 take " at “ failure rate and 

au Hr; fr f! * r “ of the glass pU ' e bei "B considered and 

an estimated total load duration for the given loading condition. The 

allowable stress for a given co mbined loading condition is the same, 
except that the allowable stress corresponding to the shortest load 
duration of the loads being considered is taken as the allowable 
stress for all loads in the combination. For example, when the 
combination of wind load (short duration) and deadweight load (long 
duration) is considered, the breakage load associated with the wind 
load is applied to the combined loading condition. 


♦The thermal stress (if any) must be determined from a different 
analysis than that described in this paper. Thermal stresses will 
exist in the glass if: 1) The glass panel is rigidly mounted t« a 

material with a different thermal expansion coefficient, 2) the 
glass panels are subjected to local heating, such as will result when 
part of the panel is shadowed from the sun, and 3) local hewing 
results from damaged solar cells attaining a back-biased condition. 
During the summer of 1979 Dave Goodwin of JPL studied the stresses 
due to local heating effects. A general conclusion of his work is 
t at for the third type of local heating effect, one should expect a 
stress in the glass of 40 lb/in 2 to 80 lb/in 2 per degree 
centigrade local temperature gradients in the glass. At this 

writing, however, no recommendations as to the temperature gradients 
to be expected are available. 
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Determine the breaking stress, Og, corresponding to each of the 
loading conditions identified in step Al. The breaking stress is a 
function of acceptable failure rate, surface area of che plate, load 
duration, and degree of temper of the glass (i.e., annealed, semi- 
tempered, tempered). The breaking stress, erg, is obtained from 
Equation 11. 

ct B " (fA><fT> CT ll (11) 


where 

a H - breaking strength of glass for 1-square-raeter samples 
subjected to a load of 1 minute total duration. This 
strength is a function of the failure rate and glass type 
and is obtained from Figure 12. 

f^ * fraction of for a plate area other than 1 square 

meter (obtained from Equation 7 or Figure 10). 

f T ■ fraction of f or a load duration other than 1 minute 

for various glass tempers (obtained from Figure 9). 

These quantities are determined as follows: 

(1) The designer must specify the acceptable failure rate. This 
should be baseJ on economic considerations; that is, the 
failure rate which produces the minimum life-cycle cost for 
a given application should be used. 

(2) From Figure 12 obtain a n> the breaking strength for the 
specified glass type for a 1-square-meter sample subjected 
to a load of 1-minute total duration corresponding to the 
failure rate chosen in step (1) above. 

(3) The fraction of for a plate area other than 1 square 
meter may be taken from Figure 10 or Equation 7, which may 
be rewritten as follows: 



where 

a ■ length of plate in meters 
b *■ width of plate in meters 

(4) Next, the designer must estimate tne total duration of the 
loading conditions delineated in step Al. Assume that the 
design life of a photovoltaic solar panel installation is 
usually 20 years, and with the exception of deadweight and 
thermal loads, the loading conditions specified in step Al 
are expected to occur only once during that period. For 
these assumptions the following ranges of total load 
duration may be recommended for the various types of loads. 


Total Load Duration 
Load for 20-Year Period 


Wind 

Earthquake 

Snow 

Deadweight 

Thermal 


1 min - 15 min 
15 sec - 5 min 
12 hr - 3 days 
20 yr 

2 yr - 7 yr 


Using Figure 9 and the above load durations, determine the 
fraction of for total load durations other than 1 

minute, depending on whether annealed, semi-tempered 
(i.e., heat strengthened) or tempered glass is being 
specified for the installation being analyzed. As 
previously mentioned, this analysis assumes that the glass 
surfaces are in good condition and have not been subjected 
to damage or "bruising 11 such as would be caused by 
mishandling, hail, impact by rocks, or sandstorms. The 
designer may interject a degree of conservatism here to 
cover these possibilities by using load durations as long 
as or longer than the higher end of the recommended^ 
range. Construct a table of load type versus fraction of 
1-minute breaking strength as shown: 


Load 

Wind 

Earthquake 

Snow 

Deadweight 

Thermal 


f 

f 

f 

f 

f 


W 

E 

S 

D 

T 


Using Equation 11 and the values of cr n , f A and f T 
obtained from steps C2, C3, and C5, respectively, the 
following table of load type versus breakage stress may be 
constructed: 


Load 


Breakage Stress 


Wind 


ct BW 


Earthquake 


17 BE 


Snow 


ct BS 


Deadweight 


a BD 


CT 


BT 


Thermal 


D. COMPARE APPLIED STRESS TO GLASS BREAKAGE STRESS 

Finally, the stresses due to the applied loads (step B5) are 
compared to the breaking stresses (step C6) to see if the assumed 
glass thickness is adequate. The glass thickness is considered 
adequate if the following combined loading criteria are satisfied. 
(Remember that for combined loading conditions the breakage stress 
associatated with the load of shortest duration considered is applied 
to the combination. ) These criteria are essentially equivalent to 
those in ANSI A58.1 -1972 (Reference 18). The criteria of other codes 
applicable to a particular installation may be substituted here. 


a BD 


1 > 

ct BL 


1 > 


a D + ( % 


ff BW ° r 


or a ) 
E 

CT BE 


1 > 


ct d 


a 3T 


1 > .75 


1 > .75 


1 > .75 


1 > .66 


i" 17 D . °L * fa W ° r °E ) ] 

L ct bw or ct be J 

[^1 

r ct d + s ° r ct e ) + ct t ] 

L ct bw or ct be J 
|\j ; ct l + (CT w or a E ) + 


a BW ° r °BE 


a 




where the "live-load" component, cr L , i s the sum of the tributory 
live load components, so that 

a L = CT S + ( an y other live loads identified) 
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SECTION VIII 


CONCLUSIONS AND RECOMMENDATIONS 


The proposed method for determining the thickness of rectangular 
glass solar collector panels is inherently very accurate. Errors in 
estimating the applied stress using this method will be canceled, 
because the same method was used to convert the extensive glass burst 
pressure data of various investigators to the glass breakage stress 
values reported herein. 

It is hoped that the proposed analytical tool will prove useful 
to the designer of photovoltaic panel installations which employ glass 
plates. An immediate need for the method exists in making life-cycle 
cost studies; for example, in the cost trade-off between decreasing 
the glass thickness, which increases the probability of failure and 
replacement costs, versus the increased first-cost and decreased 
efficiency associated with specifying thicker glass. 

In the preparation of this report the author has been sensitized 
to certain deficiencies in the proposed method. These deficiencies, 
noted below, also constitute recommendations for future work. 

(1) The actual method of determing the environmental loads is 
not within the scope of this report. It is apparent, 
however, that additional aerodynamic studies are needed to 
determine the flow field and local pressures in large 
solar array fields. Also, the thermal stresses in 
photovoltaic solar panels resulting from thermal gradients 
due to shadowing or other causes should be studied. 

(2) Total load durations to be expected over the design life 
of ^ photovoltaic panel should be studied further. 
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APPENDIX I 


SAMPLE PROBLEM 


It is desired to determine the required minimum annealed glass 
thickness for 4- foot-square photovoltaic solar panels to be installed 
in a large array of photovoltaic panels in Bangor, Maine (45° N, 

39° W) • The panels are to be designed for a 20-year life with a 2% 
failure rate at design loading conditions, Bangor, Maine was chosen 
deliberately for this sample problem, as this locality experiences 
moderately high wind, earthquake and snow loads. The outline scheme 
used below corresponds to that used in Section VII of this report, 
METHODOLOGY. 


A. DEFINE THE LOADING 

A lengthy discussion on determining the loads using ANSI A58.1 - 
1972 (Reference 18) is included here rather than in Section VII, as 
the method of determining the loads is not actually within the scope 
of this document. Furthermore, the designer may wish to use local 
building codes or other means to define the loads. 

(1) Assumed glass thickness t - .230 in. 

(2) Panel construcion and tilt angle: In addition to the 

assumed glass thickness of .230 in, it is known that the 
cells are encapsulated in a .080-in-thick layer of PVB 

( polyviny lbutyral ) and that the substrate is .032-in-thick 
aluminum. 

\ 
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(3) Define the various loading conditions in terms of uniform 
normal pressures applied to the panel. 

Wind Loads. Figure A1 of ANSI gives the basic wind speed 
(measured 30 ft above ground level) for a 25-year mean recurrence 
interval at Bangor, Maine as 


V 30 * 70 mph 
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ANSI provides additional tables to determine the pressure on signs, 
roofs, roofs over nonenclosed structures, etc., as a function of the 
geometry of the structure and the terrain (i.e., city, suburban, and 
flat, open country). These tables incorporate provisions to increase 
the dynamic pressure associated with the basic wind speed to account 
for increased pressure due to wind gusts. The tables in ANSI provide 
a reasonable means of assessing the pressure acting on photovoltaic 
panels installed on the roofs of convention buildings. Unfortunately, 
however, the aerodynamic flow in a large array of photovoltaic panels 
cannot be compared to any of the standard structures considered in 
ANSI. Pending a better understanding of the nature of this aero- 
dynamic flow, it seems desirable to derive reasonable, conservative 
values for the net local pressure acting on a photovoltaic panel by 
more direct means. This maintains a degree of visibility as to the 
assumptions included in the derivation of the wind pressure load. 

The array in which the 4-ft-square panel is mounted may be 
considered infinitely long compared to its slant height. Therefore, 
it is appropriate to employ the following graph (adapted from Hoerncr, 
Reference 19) tu determine the net normal force coefficient, C^. 



This curve shows how the net normal force coefficient on a panel 
(C^) varies with the tilt angle. The local net pressure coef- 
ficient, CpfjL, is obtained by multiplying by 2 based on the 
following argument. The center of pressure for a two-dimensional 
plate inclined to the wind is at approximately 1/3 of the chord 
measured from the windward edge of the plate. A triangular pressure 
distribution varying from zero at the trailing edge to 2 times at 
the leading edge will produce an average net normal force coefficient 
equal to C^. For this triangular pressure distribution, the center 
of pressure is at the 1/3 chord. 

These considerations result in the following relationship for 
the net local pressure P w on a photovoltaic panel which is part of a 
large array: 

P W = 2C N q 3() 
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JC. 


A 
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where 


P w - normal force per unit area due to wind to be shared by 
all layers of panel 

0 ^ = net normal force coefficient from the figure above 

^30 = dynamic pressure in lb/ft^ corresponding to the 

basic wind speed f V 30 

= .00256 (V 30 ) 2 

A degree of conservatism is introduced here due to the method of 
deriving the net local pressure coefficient. Moreover, q 3 Q 
corresponds to a height of 30 ft, and wind velocity is less at heights 
less than 30 ft. On the other hand, the present analysis is 
unconservative relative to ANSI because it does not include a gust 
factor. For the sample problem considered here 

q 3Q = .00256 ( 70 ) 2 = 12.5 lb/ft 2 

and 

p w = 2U.6K12.5) = 40.0 lb/ft 2 = .278 lb/ : .i 2 

Earthquake Loads . Equation 17 of ANSI, intended to describe the 
"lateral force on parts or portions of buildings or other structures", 
can be adapted to determine the normal pressure load on the 
photovoltaic panel as follows: 


P 


E 


P 


E 


Z 


C 


P 


w 


ZC p w 

normal force per unit area of panel surface to be 
shared by all layers of the panel 

numerical coefficient depending on the zone as 
determined from Figure A 8 (ANSI) 

Zone Z 

1 .25 

2 .50 

3 1.00 

horizontal force factor for parts or portions of 
buildings or other structures from Table 21 (ANSI) 

panel weight per unit area as calculated in A2 above 
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For this application 


Z = .50 (from Figure A8 - ANSI) 

Cp = .2 (from Table 21 - ANSI) 

w = .0289 (from step A2 above) 

So that 

P E = (.5)(.2)( .0289) = .00289 V>/in 2 

Snow Loads . ANSI provides maps with isolines of ground snow load 
in pounds per square foot for 25, 50, and 100-year mean recurrence 
intervals. For structures having no human occupants or where there is 
negligible risk to human life, the 25-year mean recurrence interval 
(Figure A7-ANSI) may be used. To find the normal pressure loads on the 
photovoltaic panels due to snow, the ground snow load is multiplied by 
a coefficient, Cg, which depends on the tilt angle, ct. The total 
pressure on the panel due to snow is 


where 


PS 


C S P GS 


P s ■ normal force per unit area of panel surface to be 
shared by all layers of the panel 

Cg = snow load coefficient 

a 

30° 

70° 

30°-70° 

P(jg s basic ground snow load for 25— year mean recurrence 
interval 


.8 

0 

.8 - 


g-30 

so 


For the present application 


P GS * 52 lb / ft2 (from Figure A7 - ANSI for Bangor, ME (45°N, 39°W)) 

(Ps^total * *5(52) ■ 26 lb/ft 2 * .181 lb/in 2 


A-4 


Deadweight Loads * The normal pressure on the panel due to 
gravity loading is obtained from the following equation: 

P D = w cos a 

For this application 

P D = .0289 cos 45° 

= .0204 lb/in 2 


Thermal Loads . Thermal loading results when temperature 
variations through the thickness, or over the surface of the glass 
induce thermal stresses. The method of defining the thermal loads for 
a particular panel design and application is not considered here. 

The total pressure loads to be shared by all of the layers of the 
panel are summarized as follows: 


Load Description 
Wind 

Earthquake 

Snow 

Deadweight 


2 

Pressure (Ib/in ) 
P W = .278 

p„ = .002 

>; ■ - 

P D - .020 


(4) Next it is necessary to determine that portion of the above 
pressure loads which the glass layer must bear. The 
cross-section of the sample problem panel is shown below: 


Glass 

PVB 

Aluminum 


t = .250 in., E - 10,000,000 lb/in 2 , v - .22 

t = .080 in., E = 300 lb/in 2 , v = .25 

t - .032 in., E » 10,000,000 lb/in 2 , v = .33 


The flexural rigidities (D) of the various layers are 


Laye r 


Et' 


12 ( 1 -/) 


Glass 13680 

PVB 

Aluminum 30 


Total 


13710 


Therefore, the fraction of the total 
glass layer is 


pressure 


loads borne by 


£l 

P 


13680 
13710 * 


.998 


The table of total 
modified to obtain 


pressure loads from 
the pressures borne 


step A3 should b^ 
by the glass layer: 


Load 

Wind 

Earthquake 

Snow 

Deadweight 


Pressure (lb/in^) 

.277 

.002 

.180 

.020 


estimate the stresses 


Now calculate the load intensity factor (LIF) f i n A 
ccrr„po„a ne atresa intensity facto/iSIF " and the str.s 
corresponding to each of the above loading ionditio's 


b = 48 in 
t = .250 in 
E = 10,000,000 lb/in 2 
v = .22 


D 


Et' 


12 (l-v Z ) 


13680 


Load 

Kind 

Earthquake 

Snow 

Deadweight 

Thermal 


.4 

LIF ■ 2 — 
Dt_ 

430 

3 

279 

31 


SIF 

from Figure 7 


a m SIF (psi) 

b Z t 


94 

2230 

1 

20 

69 

1640 

9 

210 

— 

500* 


♦Methods to calculate thermal 
value is assumed in order to 
follow* 


stress are not included here. This 
see the affect on the calculations which 


the 
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c. 


GLAb. BREAKAGE STRESS 


(1) For a failure rate of 27 . 

P f = .02 

(2) From Figure 12 we obtain the breaking strength of a 
l-square-meter sample of annealed float glass subjected to 
a load of 1-minute duration 


'11 


4200 


(3) The fraction which may be used for a 4 ft (1.22 m) square 
plate is 

f / 1 \ 1/6 

r A “ \ 1.22 x 1.22 ) 


.936 


(4) Next we estimate the total duration of the loads for a 

design life of 20 years. Remember that wind, earthquake, 
and snow loads are obtained from ANSI for a 25-year mean 
recurrence interval. For these loads, then, it remains 
only to estimate the duration of loading for one occurrence 
during the 20-year design life. The deadweight load acts 
during the entire life of the panel. The thermal load will 
obviously occur daily and may vary in intensity 
seasonally. Therefore, the total duration of the thermal 
load is probably in the range of 1/10 to 1/3 of the design 
life of the panel. The following table shows the total 
duration of the various loads, using values at the high end 
of the range recommended in Section VII, step C4. 


Total Load Duration 
Load for 20-Year Life 


Wind 

15 min 

Earthquake 

5 min 

Snow 

3 days 

Deadweight 

20 yr 

Thermal* 

7 yr 


(5) Using Figure 9 for semi-tempered glass, construct the 

following table of load type versus f T , the fraction of 
the 1-minute, l-square-meter breaking strength of annealed 
glass which may be used for other load durations. 
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Load 


f-r for Annealed Glass 


Wind 

.82 

Eartt quake 

.88 

Snow 

.62 

Deadweight 

.51 

Thermal 

.52 


(6) Finally, the 2% failure rate breakage stress for a 4 ft x 
4 ft panel of serai-tempered glass may be obtained from 
Equati n 11 for each of the load conditions considered: 

°B “ f A f T a ll (11) 

For the sample problem 

CT B = .936 f T 4200 
- 3930 f T 

The following table may now be constructed: 

ct b - 3930 f T (lb/i.i 2 ) 

Load for 2 % Failure Rate 


Wind 

3220 

Earthquake 

3460 

Snow 

2440 

Deadweight 

2000 

Thermal 

2040 


D. COMPARE APPLIED STRESS TO GLASS BREAKAGE STRESS 

The combined load criteria given in Section VII D must be 
satisfied. 



210 

2000 


CT D * 


°S 


a 


BS 


210 ♦ 1640 
2440 


.76 < 1 


g D * g W ° r a E 
a BW ° r<7 BE 


210 ♦ 2230 
3220 


- .76 < 1 
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* 


a D ^ g T 
CT BT 


210 + 500 
2040 


-, e | 210 + 1640 + 2230 | 

' [ 3220 J 

-, c I 210 ♦ 1640 + 500 | 

I 3220 J 


- .35 < 1 


= .95 < 1 


-,J a D * CT S * (CT W ° r V ] 

L °BW ° r °BE J 


.72 < 1 


.68 < 1 


The assumed thickness of .250 in. may be considered adequate for this 
application. 
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APPENDIX II 


To automate the "best-fit" Weibull analysis of brittle fracture 
data, a FORTRAN computer program was written by Robert Weaver of the 
Jet Propulsion Laboratory. A listing and sample output of this 
program follows. 


C ANALYSIS of fracture data using wiebull distribution 

r 

C CUMULATIVE PROBABILITY OF FAILURE ■ l-EXP(-( (S-SU)/SO)**MI 

c 

C WHERE ' S — MEASURED STRESS AT FAILURE - INPUT 

C SU — A STRESS BELOW WHICH NO FAILURES WILL OCCUR - OUTPUT 

C SO — A FREE PARAMETER RELATED TO SU A STRESS RANGE -OUTPUT 

c M — A POWER VALUE RELATED TO THE MATERIAL BEING SAMPLED -OUTPUT 

C A LEAST SQUARES FIT IS USED TO DETERMINE THE VALUES. THE EQUTJON USED IS 

C LOGIO ( LOG 101 1/1-CPF I ■ MLOGIO ( S-SU I -ML0G10I SO) -. 3622156 

C 

C WHICH IS ANALOGOUS TO Y« AX ♦ B WHERE* 

C 

C a * M » 8 * -MLOGIO (SO >-.3622156 (-.3622156 -LOGIOLOGIOI E J J 

C USING THE INPUTS OF STRESS AT FAILURE. S ( I > AND THE NUMBER THAT 

C FAILED AT THAT STRESS. FID AND VARYING SU TO OBTAIN A FIT CORRELATION 

C AS CLOSE TO 1 AS POSSIBLE THE PARAMETERS CAN BE DETERMINED 

C 

DIMENSION S ( 200 ) *A(200).B(200)»SU(200).SO(200)»DF(2001. 

*FS(200),TTL(12).FL(200),FLC(200) 

INTEGER F ( 200 • * TG 
DATA TTL/12*' '/ 

REAL *8 CC ( 200 ) . SCC 
FF« .3622156 
1 CALL INPT (TTL.S.F.N) 

C SORT THE INPUT FROM SMALLEST TO LARGEST 

(■;m«N-1 
K-2 

DO 15 I-l.NM 
00 10 J*K.N 

IF (S( 1 1 .LT.SIJ ) ) 60 TO 10 

TS«S( I ) 

S( I !«S< J) 

S( J 1 -TS 
TG*F 1 1 ) 

F ( I ) *F ( J ) 

F i J ) * TG 
10 CONTINUE 
X'K+1 

15 CONTINUE 
TC-0 

DSU*S ( 1 ) /200, 

DO 20 I-l.N 
20 TG-TG+FI I > 

TOTS-TG+1. 

FS ( 1 ) *F I ll/TOTS 

FL(1)»AL0G10(AL0G10(1./(1.-FS( 1 > ) I ) 

DO 30 I-2.N 

FS( II-FS < I-U+F ( 1 (/TOTS 

G-AL0G10I ALOGlOt ../ ( 1.- ( FS( I) /I . ) I )) 

FL(l)«G 
30 CONTINUE 
SSU—DSU 
DO 50 1-1.200 
SSU-SSU4DSU 

CALL FIT (SSU.SB.SA.SCC.FL.S.N) 

IF (SSU.GT .5 ( 1 ) ' GO TO 51 
A( I l-SA 
B ( I )-SB 
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CC( I ) *ScC 
SU< I ) *s$u 

SC/{ 1 ) *10*** (-( ti>B + FF ) /SA) ) 
'♦O continue 
<►1 CONTINUE 

T<>1*-ABS<CC< I i ) 

I0< 3 I 

00 50 l«2f200 
TR®1.-ABS(CC( I ) ) 

IF< TS.LT.TR I GO TO SO 
T5.TR 


50 CONTINUE 
P« At ID<) 

SGUn$U( I DK ) 

SCO-SOI I OK ) 

DO 60 I » 1 • N 

f^CU , ch7 EXP(-(<S< 1 >-SSU)/SOO)**P) 

DF ( I ) «FS < I ) -FLC ( I ) 
tO CONTINUE 

jJilf It' 6 ™) ”J-*Jo^*«s«n.Fm.Fsm.i-uNi 

•?M 

600 -™" A ' =>' 

Ml FORMA! IIHUUA./. F|, lEmiw): i "''.i 

•m,oT - & » «" 

END 

'FOR, is SUBF.SUBF 

SUBROUTINE FIT ( D . AO. A 1 ,R , Y < X ,K ) 

REAL »8 R 

REAL Y ( 200 ! «X ( 200 ) , XX ( 200 ) 

sx-o. 

SY*0. 

sxs-o. 

SYS-0. 

SXY»0. 


DU 9 J-l ,X 

9 XX( I )"ALOG10(X( I)-D) 
DO 10 I ■ 1 . X 
SX -SX+XX1 I ) 

SXS»SX6 + XX( I ) * XX l I ) 
SY-SY + YI I ) 

SYS-SYS+YI I ) *Y ( I | 

SXY -SXY +XX ( I )*Y( I ) 
10 CONT INUE 

Z* K*5XY-SX»SY 
*• K*SXS-SX*SX 
Al-Z/kf 


A0-ISY-A1*5X)/X 

R.- Z/SQR T (W*|X*SYS-SY*SY) ) 

RE ’URN 

END 

'FOR. IS SUB I .SUB I 

SUBROUTINE INPT (T.S.F.n) 
DIMENSION T ( 1 2 ) . s ( 200 ) 
INTEGER F I 200 ) 

READ (5, 100. END-99) T 
100 FORMAT ( 12A6 ) 


>3.iUl) N.IS(I) ,FI I I . t, 
101 FORMAT ( I5/(6IF6.2.14J>) 
RFTURN 
99 STOP 
END 


• XOT 


BOWLES AND SUGARMAN - .197 PLATE GLASS 
10 

*301 1 4.675 2 5.039 5 5,401 

7.135 4 7.470 1 


3 


5.754 


3 6*106 4 6.452 3 6.795 


4 
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SAMI’ 1 .1'- OUTPUT 


BOWLES AND SUGARNaN - .197 PLATE GLASS 
INPUT DATA 


SAMPLE SIZE • 31. 


STRESS 

HJ FAIL 

PROS 

4.3010 

1 

• 0323 

5.4010 

3 

.3346 

6.4520 

3 

.6774 

7.4700 

1 

.9677 


STRESS 

NO MIL 

PROS 

4.6750 

2 

.0968 

5.7540 

3 

.4516 

6.7950 

4 

.6065 


stress 

NO FAIL 

PROB 

5.0390 

3 

.2581 

6.1060 

4 

.5606 

7.1350 

4 

.9355 


BOWLES AND SUGARMAN - .197 PLATE GLASS 
FIT RESULTS 


su 

SLOPE 

so 

•00000 

7.74879 

6.25482 

• 021 30 

7.71961 

6.23304 

•04301 

7.69039 

6.21127 

• 06451 

7.66120 

6 . 18949 

.08602 

7.63199 

6.16770 

. 10752 

7.60278 

6. 14502 

• 12903 

7.57356 

6.12413 

• 15053 

7.54435 

6.10234 

. 17204 

7.51514 

6.08055 

. 19354 

7.48592 

6.05876 

•21505 

7.45666 

6.03696 

•23655 

7.42742 

6.01516 

.25806 

7.39818 

5.99336 

. 27956 

7,36892 

5.97156 

.30107 

7.33966 

6.94976 

.32257 

7.31041 

5.92795 

. 34408 

7.28113 

5.90614 

. 36558 

7.25184 

5.88433 

.38709 

7.22257 

5.86251 

•40859 

7.19326 

5.84070 

.43010 

7.16398 

5.81888 

•45160 

7.13466 

6. 79706 

•47311 

7.10535 

6. 77523 

.49461 

7.07602 

5.75341 

.5161? 

7.04669 

5 • 73158 

.53762 

7.01737 

5. 70974 

.55913 

6.98802 

5.68791 

.58063 

6.95867 

6.66607 

• 6 02 1* 

6.92932 

5.64423 

• 62 364 

6.69996 

5.62239 

.64515 

6.67059 

5.60054 

.66665 

6.84121 

5.57869 

•66816 

6.61182 

5.55684 

. 70966 

6.78243 

5.53499 

•73117 

6.75302 

5.51313 

• 75267 

6 . 72 36 1 

5.49127 

•77418 

6 . 694 1 9 

5.46941 

.79568 

6.66476 

5.44754 

•81719 

6.63531 

5.42567 

.83869 

6.60588 

5.40380 

.86020 

6.57641 

5.38192 

.88170 

6.54695 

5.36004 

.90321 

6.51747 

5.33816 

.92471 

6.48799 

5.31628 

.94622 

6.45850 

5.29439 

.96772 

6.42900 

5.27249 

.98923 

6.39949 

5.23060 


CORRELATION COEFFICIENT 

.9826088960+00 

.9826511821+00 

.9826919805+00 

. 9627 3A 5633+00 

.9827773651+00 

.9826199653+00 

.9826630257+00 

.9829070851+00 

.9829520528+00 

.9829969426+00 

.9830396806+00 

.9830854277+00 

.9831305015+00 

.9831757920+00 

.9832218114+00 

.9632686015+00 

.9833148776+00 

.9633604366+00 

.9834086231+00 

.9834549674+00 

.9835040349+00 

.9835511962+00 

.9835996009+00 

.9836471679+00 

.9834961522+00 

.9837461611+00 

.9837956245+00 

.9838452970+00 

.9838958279+00 

.9839474415+00 

.9839978082+00 

.9840492998+00 

.9841008224+00 

.9841533728+00 

.9842056210+00 

.9842582531+00 

.9843119056+00 

.9843648696+00 

.9844184211+00 

.9844735496+00 

.9645274450+00 

.9845627364+00 

.9846377137+00 

.9846939970+00 

.9847503351+00 

.9848073403+00 

.9846643536+00 
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A 


1.01073 

6.36997 

5.22870 

•9849218346+00 

1.03224 

6.34044 

5.20680 

.9849600075+00 

1.05374 

6.31089 

5.18489 

.9850378863+00 

1.07525 

6.28134 

5.16298 

.9850965493+00 

1.09675 

6.25179 

5. 14106 

.9851567372+Ot) 

1.11826 

6.22221 

5.11915 

.9852160854+00 

1.13976 

6.19262 

5.09723 

.9852756053+00 

1. 16127 

6.16303 

5.07530 

.9853359298+00 

1.16277 

6.13343 

5.05337 

.9853978754+00 

1.20428 

6.10381 

5.03144 

.9854584757+00 

1.22576 

6.07418 

5.00950 

.9855204595+00 

1.24729 

6.04454 

4.98756 

.9855830103+00 

1.26879 

6.01489 

4.96561 

.9856462789+00 

1.29030 

5.98523 

4.94366 

.9857091161+00 

1.31180 

5.95555 

4.92171 

.9857731124+00 

1.33331 

5.92587 

4.89975 

.9858378088+00 

1.35481 

5.89616 

4.87779 

.9859022688+00 

1.37632 

5.86645 

4.85582 

.9859674504+00 

1.39782 

5.83672 

4.83385 

.9860334153+00 

1.41933 

5.80698 

4.81187 

.9860992245+00 

1.44083 

5.77723 

4.78989 

.9861666600+00 

1.46234 

5.74746 

4.76791 

.9862332157+00 

1.48384 

5.71768 

4.74591 

.9863010625+00 

1.50535 

5.68788 

4.72392 

.9863692303+00 

1.52685 

5.65807 

4.70192 

.9864379772+00 

1.54836 

5.62824 

4.67991 

.9865074156+00 

1.56936 

5.59840 

4.65790 

.9865774662+00 

1.593 37 

5.56855 

4.63589 

.9866477199+00 

1.61287 

5.53867 

4.61387 

.9867179734+00 

1.63438 

5.50879 

4.59184 

.9867906709+00 

1.65588 

5.47889 

4.56981 

.9868625541+00 

1.67739 

5.44897 

4.54777 

.9869354033+00 

1.69889 

5.41903 

4.52572 

.9870086336+00 

1. 72040 

5.38907 

4.50368 

.9870817242+00 

1.74190 

5.35910 

4.48162 

.9871562065+00 

1.76341 

5.32912 

4.45956 

.9872313189+00 

1.78491 

5.29911 

4.43749 

.9873071705+00 

1.80642 

5.26909 

4.41542 

.9873833679+00 

1.82792 

5.23904 

4.39334 

.9874593557+00 

1.84943 

5.20898 

4.37125 

.9875371575+00 

1.87093 

5.17889 

4.34916 

.9876144769+00 

1.89244 

5.14879 

4.32706 

.9876934116+00 

1.91394 

5.11867 

4.30496 

.9877725639+00 

1.93545 

5.08853 

4.28285 

.9878524562+00 

1.95695 

5.05837 

4.26073 

.9879328459+00 

1.97846 

5.02818 

4.23860 

.9880136680+00 

1.99996 

4.99797 

4.21647 

.9880955034+00 

2.02147 

4.96774 

4.19433 

.9881776085+00 

2.04297 

4.93749 

4.17218 

.9882606918+00 

2.06448 

4.90721 

4.15002 

.9883443251+00 

2.08598 

4,87691 

4,12786 

.9884284999+00 

2. 10749 

4.84659 

4.10569 

.9885136722+00 

2. 12899 

4.81624 

4.08351 

.9885989108+00 

2. 15050 

4.78587 

4.06132 

.9886853658+00 

2. 17200 

4.75547 

4.03913 

.9887724110+00 

2. 19351 

4.72504 

4.01692 

.9888598558+00 

2.21501 

4.69458 

3.99471 

.9889478307+00 

2.23652 

.66410 

3.97249 

.9890370207+00 

2.25P02 

4.63359 

3.95026 

.9891260914+00 

2.27953 

4.60305 

3.92802 

.9892169553+00 

2.30103 

4.57248 

3.90577 

.9893078882+00 

2.32254 

4.54188 

3.88351 

.9893997019+00 

2.34404 

4.51125 

3.86125 

.9894919767+00 

2.36555 

4.48059 

3.83897 

.9895849262+00 

2.38705 

4.H4989 

3.81668 

.9896787159+00 

2.40856 

4.41916 

3.79438 

.9897730695+00 

2.43006 

4.38840 

3.77207 

.9898684225+00 

2.45157 

4.35760 

3.74975 

.9899642513+00 

2.47307 

4.32677 

3.72742 

.9900608871+00 


2*49458 

4.29590 

3.70508 

. 9901580357+00 

2.51608 

4 • 26499 

3.68273 

. 9902559371+00 

2.53759 

4 . 234 C 4 

3.66036 

. 9903542879+00 

2.55909 

4.20305 

3.63799 

. 9904535904+00 

2.50060 

4.17202 

3.61560 

. 9905534991+00 

2.60210 

4.14095 

3.59320 

. 9906541308+00 

2.62361 

4.10984 

3.57079 

. 9907553136+00 

2.64511 

4.07868 

3.54836 

. 9908569458+00 

2.66662 

4.04747 

3.52592 

. 9909593801+00 

2.68812 

4.01622 

3.50347 

. 9910624927+00 

2-70963 

3.98492 

3.48100 

. 9911662636+00 

2.73113 

3,95357 

3.45852 

. 9912702959+00 

2 . 7526 * 

3.92217 

3.43602 

. 9913753551+00 

2.77414 

3.89072 

3.41351 

. 991480581 7+00 

2.79565 

3.85921 

3.39099 

. 9915865252+00 

2.81715 

3,82765 

3.36844 

. 9916929809+00 

2.63866 

3.79603 

3.34589 

. 9917997859+00 

2.86016 

3.76435 

3.32331 

. 9919069037+00 

2.88167 

3.73261 

3.30072 

. 9920148875+00 

2.90517 

3.70081 

3.27812 

. 9921229264+00 

2.92468 

3.66894 

3.25549 

* 992231 38 93 +OC 

2.94618 

3.63701 

3.23285 

. 9923401197+00 

2.96769 

3.60500 

3.21018 

. 9924490486+00 

2.98919 

3.57293 

3.18750 

. 9925582410+00 

3.01070 

3.54078 

3.16480 

. 9926674987+00 

3.03220 

3.50855 

3.14208 

. 9927767582+00 

3.05371 

3.47624 

3.11934 

. 9928859936+00 

3.07521 

3.44385 

3.09658 

, 9929952352+00 

3.09672 

3.41138 

3.07380 

. 9931042411+00 

3 . 1182 2 

3.37682 

3.05099 

• 9932 1 302 35 + 00 

3 . 13973 

3.34616 

3.02017 

. 9933213619+00 

3 . 16123 

3.31341 

3 • 005 32 

. 9934291526+00 

3 . 18274 

3.28057 

2.98244 

. 9935364722+00 

3 . 2042 ^ 

3.24761 

2.95954 

. 9936 ^ 30695+00 

3.22575 

3.21456 

2.93662 

. 9937486504+00 

3.24725 

3.18139 

2.91366 

. 9938531873+00 

3.26876 

3.14810 

2.89069 

. 9939564455+00 

3.29026 

3.11469 

2.86768 

. 9940582734+00 

3.31377 

3.08116 

2.84465 

. 9941583949+00 

3.33327 

3.04750 

2.82158 

. 9942565980+00 

3.35478 

3.01369 

2.79649 

. 9943525730+00 

3. 37628 

2.97975 

2.77536 

. 9944460561+00 

3.39779 

2.94565 

2.75221 

. 9945367513+00 

3.41929 

2.91140 

2.72902 

. 9946241388+00 

3.44080 

2.87697 

2.70579 

. 9947079968+00 

3.46230 

2.84238 

2.68253 

. 9947877330+00 

3.48581 

2 . 8 U 760 

2.65924 

. 9948629166+00 

3.50531 

2.77263 

2.63591 

. 9949329690+00 

3.52682 

2.73745 

2.61254 

. 9949972191+00 

3.54832 

2.70207 

2.58913 

. 9950550084+00 

3.56983 

2.66645 

2.56567 

. 9951055172+00 

3.59133 

2.63059 

2.54218 

. 9951479208+00 

3.61284 

2.59448 

2.51864 

. 9951811365+00 

3.63434 

2.55810 

2.49506 

. 9952041033+00 

3.65585 

2.52143 

2.47144 

. 9952154964+00 

3.67735 

2.48445 

2.44776 

. 9962139284+00 

3.69886 

2.44715 

2.42404 

. 9951976830+00 

3.72036 

2.40949 

2.40027 

. 9951650182+00 

3.74187 

2.37146 

2.37644 

. 9951136272+00 

3.76337 

2.33302 

2.35257 

. 9950411247+00 

3.78485 

2.29415 

2.32864 

. 9949446769+00 

3.80638 

2.25480 

2.30465 

. 9948209543+00 

3.82789 

2.21494 

2.28061 

. 9946661003+00 

3*84939 

2.17453 

2.25652 

. 9944756377+00 

3.07090 

2.13350 

2.23237 

. 9942442641+00 

3.69240 

2.09160 

2.20816 

. 9939656998+00 

3.91391 

2.04936 

2.18390 

. 9936325087+00 

3.93541 

2.00611 

2.15959 

. 9932356766+00 

3.95692 

1.96194 

2.13524 

. 9927643490+00 


B -5 


3.97642 
3.99993 
4.02143 
4.04294 
4.06444 
4.08595 
4. 10745 
4. 12896 
4. 15046 
4.17197 
4. 19547 
•..21498 
4.23648 
4.25799 
4.27949 


1.91675 
1 .87042 
1.82277 
1.77361 
1.72271 
1.66977 
1.61439 

1.55605 

1.49407 

1.42744 

1.35470 

1.27352 

1.17981 

1.06521 

.90646 


2.11084 

2.06641 

2.06196 

2.03751 

2.01309 

1.98873 

1.96450 

1.94049 

1.91683 

1.89376 

1.87166 

1.85127 

1.83408 

1.82386 

1.83421 


.9922052438+00 

.9915420016+00 

.9907541233+00 

.9898157370+00 

.9886936416+00 

.9873443398+00 

.9857098778+00 

.9837108586+00 

.9812353388+00 

.9781191662+00 

.9741091457+00 

.9687873500+00 

•9613913861+00 

.9502896941+00 

.9307063476+00 


BOWLES AND SUGARMAN - .197 PLATE GLASS 


USING CORRELATION COEFFICIENT OF .9952154964+00 
su ■ 3.656 SO « 2.471 SLOPE « 2.5214 


STRESS NO FAIL 
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